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1. Introduction 

Strongly interacting matter undergoes a deconfining transition at some temperature wiiich is 
triggered by large increase in the number of degrees of freedom as well as melting of hadronic 
states. Another important feature of the deconfined phase is the screening of color charges. It 
has been argued that color screening will lead to quarkonium dissociation above the deconfine- 
ment temperature which can be used as a signature of quark gluon plasma formation in heavy ion 
collisions [[I]]. Melting of quarkonium states has been studied in potential models with screened 
potentials (see e.g. ^ and references therein). Alternatively, one may try to reconstruct quarko- 
nium spectral functions from Euclidean correlators (see e.g. and references therein). It turns 
out, however, that quarkonium melting does not affect the Euclidean meson correlators and these 
correlators are almost temperature independent. This is contrary to what happens in the light quark 
sector, where both meson correlators and spectral functions show significant temperature depen- 
dence in the deconfined phase [||]. Therefore potential models are useful tools to study in-medium 
quarkonium properties. But to establish their applicability a better understanding of color screening 
is needed. 

On the lattice color screening is usually studied in terms of the Polyakov loop correlator related 
to the free energy of static quark anti-quark pair [^]. This correlator shows significant temperature 
dependence across the transition and in the deconfined phase the free energy of static quark anti- 
quark pair shows large temperature dependence even for very small separations between the static 
quark and anti-quark, much smaller than the inverse temperature. In perturbative picture this can 
be understood due to the fact that in the deconfined phase not only singlet quark anti-quark {QQ) 
states contribute to the free energy but also colored states with QQ in the adjoint representation. 
This observation is also supported by lattice calculations of the correlation function of two temporal 
Wilson lines in Coulomb gauge, which in perturbation theory corresponds to the so-called singlet 
free energy [^]. The singlet free energy is temperature independent at short distances and coin- 
cides with the zero temperature potential as expected. At large distances, however, it also shows 
significant temperature dependence. 

The problem of defining color singlet and adjoint QQ states on the lattice has been considered 
in Ref. [^. It has been found that the conventional definition of singlet and adjoint states have 
problems. Here we report on our study ^ of static meson correlators in 4 dimensional SU (2) 
gauge theory at finite temperature and show how the problem observed in Ref. can be resolved 
in the limit of small distances and/or high temperatures. 

2. Static meson correlators 

On the lattice correlators of the static meson operators in color singlet and adjoint states at 
t = \/T have the form [0]: 

G,{r,T) = ^^{i:vL\x)U{x,y,0)L{y)U\x,yMT)), r=\x-yl (2.1) 

Ga{r,T) = ^^^{TvL\x)TvL{y)) - -^—lTvL\x)U{x,y-Q)L{y)U\x,y, 1/r)), (2.2) 
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where N is the number of colors {N = 2m our numerical calculations). Here L{x) is the temporal 
Wilson line, which on the lattice is simply L = n^io^ Uq{x,x) with Uq{x,x) being the temporal 
links. The correlators depend on the choice of the spatial transporters U{x,y;t). In the special 
gauge, where U {x,y;t) = 1 the above correlators give standard definition of the singlet and adjoint 
free energies 

,-F,{nr)/r^l^TrLt(x)L(,)), ^ (TrLt(.)TrL(y)) _ (TrLt(.)L(y)) 

The singlet and triplet free energies can be calculated at high temperature in leading order HTL 
approximation [^. At this order F\ and Fa are gauge independent or in other words do not depend 
on the choice of the parallel transporters U{x,y;t). At small distances the singlet free energy is 
temperature independent and coincides with the zero temperature potential, while the adjoint free 
energy depends on the temperature [Q. 

The physical free energy of a static quark anti-quark pair is given by the thermal average of 
the singlet and adjoint free energies and is explicitly gauge independent 

e-'^'-'^" = ^--"'^^'^^/^ + ^--""^^'^^/^ = ^(TrL(.)TrL(,)) = i,G(r,r). (2.4) 
Using the ti^ansfer matrix one can show that in the confined phase 

oo oo 

Gi(r,r) = £c„(r)^'-^"('-'^)/^ G{rJ) = £ ^-^"('•-n/r^ (2.5) 

where are the energy levels of static quark and anti-quark pair. The coefficients c„ (r) depend 
on the choice of the gauge transporters. The color averaged correlator G{r,T) does not contain 
c„. The lowest energy level is the usual static quark anti-quark potential, while the higher energy 
levels correspond to hybrid potentials [[lO|]. If ci = 1 the dominant contribution to Ga would be 
the first excited state E2, i.e. the lowest hybrid potential which at short distances is related to the 
adjoint potential. In this sense Ga is related to static mesons with QQ in adjoint state. Numerical 
calculations show, however, that c\ (r) / 1 and depends on the separation r. Thus Ga also receives 
contribution from Ei [^]. The lattice data seem to suggest that ci approaches unity at short distances 
in accord with expectations based on perturbation theory, where ci = 1 up to ^(g^) corrections 



[II]. Therefore at short distances, r <^\/T the color singlet and color averaged free energy are 
related F(r, T) = Fi{r, T) + T\n{N'^ — 1)- This relation is indeed confirmed by lattice calculations 



[12]. 



3. Numerical results 

We have calculated correlation functions of static mesons Gi{r,T) and G{r,T) both in the 
confined and deconfined phase of SU (2) gauge theory. The details of the calculations are presented 



in Ref. [jSp. We have studied the color singlet and averaged correlators given by Eqs. ( [2.1| ) and ([2.4|). 
The spatial links entering the transporter U {x,y;0) were smeared using APE smearing, which has 
been applied iteratively. The weight of the staple in the APE smeared link was 0.12. For j8 = 2.5 
we use spatial links with 10 steps of APE smearing and unsmeared spatial links. For j3 = 2.7 we 
used unsmeared spatial links as well as spatial links with 10 steps and 20 steps of APE smearing. 
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Figure 1: The color averaged free energy below 
deconfinement temperature at j3 = 2.5 calculated 
on 32^ X Nx lattices. Also shown is the T = po- 
tential. 
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Figure 2: The color singlet free energy below de- 
confinement temperature at ji — 2.5 calculated on 
32^ X Nr lattices. It is also shown with the contri- 
bution rinci subtracted (inset). 



3.1 Color averaged correlator in the confined phase 

The color averaged correlator has been calculated in the confined phase in the temperature 
interval 0.327^ - 0.957^ for j8 = 2.5 and 0.497^ - 0.987^ for jS = 2.7. The numerical results for 
the color averaged free energy for j8 = 2.5 are shown in Figure To ehminate the trivial tem- 
perature dependence due to the color trace normalization in Figure [I] we show the subtracted free 
energy F'{r, T) =F{r,T) — T\nA together with zero temperature potential. The color averaged free 
energy does not show any temperature dependence up to temperatures of about 0.767^. Since the 
temperature dependence for T < 0.16Tc is relatively small we attempted to fit the color averaged 
correlator with the 1-exponential form G{r,T) = c"{r)exp{—Ei{r) /T). The ground state energy 
E\ (r) extracted from this fit agrees well with the zero temperature potential calculated in Ref. [p^], 
while the coefficients Cj(r) are close to one as expected. 

3.2 Color singlet correlators in the confined phase 

The color singlet correlators have been calculated using different levels of APE smearing in 
the spatial gauge connection. We have found that when no smearing is used the color singlet free 
energy, —TlnGi{r,T) shows a small but visible temperature dependence. In particular F\{r,T) 
is larger than the T = potential for intermediate distances 0.5 < ry/o < 2. The temperature 
dependence of the singlet free energy is significantly reduced when APE smearing is applied. In 
Figure]^ we show the color singlet free energy for jS = 2.5 and 10 APE smearings. As one can see 
from the figure the color singlet free energy shows much smaller temperature dependence as we 
get closer to the deconfinement temperature. 

To understand the temperature dependence of the color singlet correlator we use 1-exponential 
fit G\ {r,T) = ci (r) exp(— £"1 {r)/T). In all cases considered the values of Ei (r) extracted from fits 



are in good agreement with the calculation of the zero temperature potential in Ref. [10]. The 
value of the prefactor c\ (r) is shown in Figure ^ When no APE smearing is used the value of ci (r) 
strongly depends on the separation r. At small distances it shows a tendency of approaching unity as 
one would expect in perturbation theory and decreases with increasing distance r. At large distance 
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Figure 3: The pre-exponential factor of the color singlet correlators as function of distance r for j3 = 2.5 
(left) and j3 = 2.7 (right). Shown are results for unsmeared spatial links and 10 and 20 steps of APE smearing. 



its value is around 0.3 — 0.5. Similar results for c\ have been obtained in study of SU (2) gauge 
theory in 3 dimensions [^]. When APE smearing is applied the r-dependence of the amplitude 
c\ is largely reduced and its value is close to unity both for j8 = 2.5 and j8 = 2.7. For j8 = 2.7 
we also see that increasing the number of smearing steps from 10 to 20 reduces the deviation of 
c\ from unity. In any case at sufficiently short distances ci is very close to one as expected in 
perturbation theory. Thus, almost the entire temperature dependence of the singlet free energy at 
distances 0.5 < r^/o < 2 is due to the deviation of c\ from unity and can be largely reduced by 
applying APE smearing to the Unks in the spatial gauge connections. To further demonstrate this 
point in the inset of Figure § we show the results for Fi (r, T) + Tlnci (r). 

3.3 Color singlet free energy in the deconfined phase 

It turns out that the singlet free energy, F\ {r,T) = —TlnGi (r, T), calculated from cyclic Wilson 
loops shares the same qualitative features as the singlet free energy calculated in Coulomb gauge 
[jl3|]. At short distances it is temperature independent and coincides with the zero temperature 
potential. At large distances it approaches a constant Foo{T), which is the free energy of two 
isolated static quarks at infinite separation. 

At leading order F\{r,T) — Foo{T) is of Yukawa form, therefore in Figure ^ we show our 
numerical results in terms of the screening function S{r,T) = r- {F\{r,T) — Foo{T)) at different 
temperatures. At short distances (rT < 0.5) the singlet free energy does not depend on the smearing 
level. Furthermore, it is very close to the free energy calculated in Coulomb gauge. We expect that 
at large distances the screening function S{r,T) will show an exponential decay determined by a 
temperature dependent screening mass m\{T), which is equal to the leading order Debye mass 
up to the non-perturbative corrections: mi = mo + ^{g^)- From Fig. ^ we can see that indeed 
S{r, T) behaves exponentially with screening mass proportional to the temperature. Fitting the large 
distance behavior of the screening function by an exponential form exp(— mi(r)r) we determine 
the screening mass m\{T). In Fig. ^ we also show the color singlet screening masses extracted 
from the fits and compare them to the results obtained in Coulomb gauge in Ref. [ pj| ] as well as to 
the leading order Debye mass calculated using 2-loop gauge coupling g{pi = 2nT) in M5-scheme. 
As we see from the figure the screening masses are smaller than those calculated in Coulomb gauge 
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Figure 4: The screening function at different tem- 
peratures as function of rT and different values of 
j3. In the inset the screening masses m\ extracted 
from singlet free energies are shown. The line 
shows the leading order results for the Debye mass. 
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Figure 5: The triplet free energy at different tem- 
peratures calculated at j3 = 2.5. The filled sym- 
bols correspond to calculations in Coulomb gauge. 
Also shown is the first hybrid potential calculated 
inRef. 



and agree well with the leading order perturbative prediction. 
3.4 Color triplet free energy 



We have calculated the color triplet correlator defined by Eq. ( |2.2| ) for different temperatures 
below and above the transition temperature. Below the deconfinement temperature we observe a 
moderate T-dependence of the triplet correlator. We also find that the corresponding free energy 



-T\nGT,{r,T) is smaller than the first hybrid potential calculated in Ref. [10], but larger than the 



triplet free energy in Coulomb gauge [13] 



Let us assume that only two states contribute to the Eqs. (p^. Then from Eq. ( p^ ) it follows 

that 

F3(r,r)=£2(r)-rin(l-C2(r) + l(l-ci(r))^.'^W/^), (3.1) 

withA£'(r) =E2{r)—E\{r) andc2 <C 1. Wehave subtracted the correction rin(l + ^(l — ci)e^^/^) 
from the triplet free energy assuming that E\ {r) is given by the ground state potential and Ej {r) is 
given by the first hybrid potential as calculated in Ref. [[I0[]. The numerical results are summarized 
in Fig. H which shows that after this correction is accounted for in the confined phase the triplet free 
energy at low temperatures agrees reasonably well with the first hybrid potential. As temperature 
increases more excited states contribute. In particular, at 0J6Tc the value of the triplet free energy 



can be accounted for by including the next hybrid state [[rO|]. 
4. Conclusions 

We have studied singlet and triplet static quark anti-quark correlators in finite temperature 
SU (2) theory expressed in terms of Polyakov loop correlators and cyclic Wilson loops. In leading 
order and probably next-to-leading order of perturbation theory the static correlators defined by 
Eq. ( |2. 1[ ) and (2.2) project onto singlet and triplet states respectively, however, this separation 



does not hold in general case. Due to interactions with ultrasoft fields there will be a mixing of 
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singlet and triplet states which is proportional to a^{\/r) and {rjx)'^, with ji being the ultrasoft 
scale JTl]]. In our case the ultrasoft scale can be the binding energy, a^/r, Kqqd or . Therefore 
it is expected that mixing is quite small at sufficiently small distances. We determined the mixing 
between singlet and triplet states in terms of the overlap factor c\(r). If the overlap factor is unity 
there is no mixing. Our lattice calculations show that c\ indeed approaches one at small distances. 
Therefore the contribution of singlet state to G-i(r^T) appears to be small at temperatures close 
to deconfinement temperature. This contribution is also controlled by the non-perturbative gap 
between the singlet and triplet states, i.e. the gap between the static potential and the first hybrid 
potential. 

Our analysis shows that at short distances rJ < 1 the singlet correlator is almost temperature 
independent, while the triplet correlator is largely affected by the deconfinement. The temperature 
dependence of the triplet correlators indicate the melting of the non-perturbative gap between the 
singlet and the triplet states above deconfinement, which turns out to be consistent with perturbative 



expectations. This finding is important for application of thermal pNRQCD discussed in Ref. [ |14{ ] 
to realistic quarkonia and temperatures not very far from the deconfinement temperature. 
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